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Abstract 

The genus Dolichopoda (Orthoptera; Rhaphidopohoridae) is present in Italy with 9 species distributed 
from northwestern Italy (Piedmont and Liguria) to the southernmost Apennines (Calabria), occurring 
also in the Tyrrhenian coastal areas and in Sardinia. Three morphologically very close taxa have been de- 
scribed in Piedmont and Liguria, i.e., D. ligustica ligustica, D. ligustica septentrionalis and D. azami azami. 

To investigate the delimitation of the northwestern species of Dolichopoda, we performed both mor- 
phological and molecular analyses. Morphological analysis was carried out by considering diagnostic char- 
acters generally used to distinguish different taxa, as the shape of epiphallus in males and the subgenital 
plate in females. Molecular analysis was performed by sequencing three mitochondrial genes, 12S rRNA, 
16S rRNA, partially sequenced and the entire gene of COL 
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Results from both morphological and molecular analyses highlighted a very homogeneous group of 
populations, although genetically structured. Three haplogroups geographically distributed could be dis- 
tinguished and based on these results we suggest a new taxonomic arrangement. All populations, due to 
the priority of description, should be assigned to D. azami azami Saulcy, 1893 and to preserve the names 
ligustica and septentrionalis, corresponding to different genetic haplogroups, we assign them to D. azami 
ligustica Baccetti & Capra, 1959, stat. n. and to D. azami septentrionalis Baccetti & Capra, 1959, stat. n. 

Keywords 

Westernmost Dolichopoda species, barcoding, systematics 



Introduction 

Dolichopoda is a Mediterranean genus, consisting of more than fifty species distributed 
throughout the North Mediterranean regions from the Pyrenees to Transcaucasia and 
north Iran (Di Russo and Rampini 2014). The highest species diversity is present in in- 
sular and peninsular Greece where most of the species have been described (Di Russo et 
al. 2007). However, many species are present also in Italy. In particular, nine species of 
Dolichopoda occur in Italy, where they range from the Maritime Alps to the southern tip 
of the Italian peninsula, excluding the oriental Alps, the Liguria Apennines, the Apuane 
Alps and the southern Apulia. Most species of this genus are strictly dependent upon 
caves. However, especially in the northern part of the range, several populations inhabit 
cave-like habitat such as rock-crevices and ravines and individuals can be observed, 
during the night, outside in moist or mesic woods. In peninsular Italy, Dolichopoda 
populations often live in cellars, catacombs, aqueducts, Etruscan tombs and other simi- 
lar man-made hypogean environments. Population sizes can be small and constant over 
long periods, at least in natural caves (Carchini et al. 1983; Sbordoni et al. 1987). 

Based on morphology three subgenera have been described in Italy: Dolichopoda 
{Dolichopoda), D. {Chopardina) and D. {Capraiacris) (Baccetti 1975; Baccetti and 
Capra 1959, 1970). One of the main morphological differences between these sub- 
genera is presence or absence of spinulation on the legs. The subgenus Dolichopoda 
includes the highest number of species distributed throughout the range of the genus, 
except for some coastal areas and it is characterized by the presence of spines on the 
anterior tibiae. Species belonging to Chopardina sub-genus are mostly found in insular 
and peninsular Tuscany and Sardinia in Italy and in Corsica, with two other species 
found in Greece (Macedonia and Etolia Akarnania). Chopardina species are character- 
ized by the presence of several spines also on the hind femurs. The subgenus Capraiacris 
includes only two species, restricted to the Giglio Island and Monte Argentario in the 
Tuscan archipelago. Its two species {D. aegilion and D. baccettii) are distinct from all 
the other species because of the lack of spines on the anterior tibiae and on the hind fe- 
mur. However, studies based on molecular characters did not support the subdivision 
in sub-genera, showing that the subgenus Chopardina is polyphyletic and the subgenus 
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Dolichopoda is paraphyletic (Allegrucci et al. 2005). Since the main morphological 
character defining the subgenera is the presence/absence of spines on the appendages, 
Allegrucci et al. (2005) concluded that such character is an adaptive character subject 
to homoplasy. Based on these observations we decided to not consider the subdivision 
in sub-genera of the species belonging to genus Dolichopoda. 

Based on molecular phylogenetic reconstruction (Allegrucci et al. 2005; 2009; 
2011) the current distribution of Dolichopoda species has been essentially shaped by 
the palaeogeographic and climatic events occurred in the Mediterranean region, start- 
ing from Late Miocene. However, also adaptation to cave life seems to have played an 
important role in this process. In particular, the end of Messinian salinity crisis would 
mark the geographic separation of the epigean forest populations and the beginning of 
dispersal toward the west Mediterranean region. The migration and /or the isolation 
of the different insular populations might have also been favored by marine regressions 
and transgressions, during the Plio-Pleistocene. This era has been also characterized by 
climatic shifts during alternating glacial and interglacial periods. The sylvicolous ances- 
tors of Dolichopoda might have used caves as refugia during the unfavorable climatic 
conditions, beginning their adaptation to subterranean habitat. Therefore, the current 
distribution of Dolichopoda can be explained by a combination of both vicariance and 
dispersal events, with many processes occurring in ancestral epigean populations be- 
fore the invasion of the subterranean environment. 

The most western Italian species are represented by D. ligustica, and D. azami. 
Two subspecies are attributed to D. ligustica, i.e., D. ligustica ligustica and D. I sep- 
tentrionalis. The first one is widespread in the Southern Piedmont (Maritime Alps) 
and in the north-western Liguria while the range of D. I septentrionalis is northern- 
most and limited to a small area in Val di Lanzo valley, close to Turin, and in some 
easternmost caves, close to Bergamo. D. azami Saulcy, 1893 is widespread in the 
south eastern France and observed also in Piedmont in a cave in Grana Valley (Bac- 
cetti and Capra 1959). The systematic position of D. azami has long been discussed 
(Finot 1890; Azam 1893; Chopard 1922, 1951) until Baccetti (1966) reclassified 
the systematics of the western Dolichopoda species, based on the morphology of the 
epiphallus and of the tergum 10, the most important characters for the systematics of 
the genus Dolichopoda. 

D. ligustica and D. azami have very similar morphologies with little variation in 
the characters on which the systematics are based. To better understand the differ- 
entiation between these two species and the limits of the ranges occupied by each of 
them, we carried out a study on several populations of D. ligustica and D. azami, both 
at morphological and molecular levels. The main aims of this study are: i) to infer the 
phylogenetic relationships of these two species, using molecular markers; ii) to com- 
pare their differentiation with the differentiation of all other species belonging to genus 
Dolichopoda previously analyzed (Allegrucci et al. 2005, 2009, 201 1); iii) to revise the 
taxonomic arrangement of species and subspecies based on both morphological and 
molecular characters. 
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Figure I . Sampling sites of Dolichopoda populations considered in this study (see also Table 1). Triangles 
indicate the type locality of the three considered taxa (CHA: D. azami, TOR: D. I. ligustica, PUG: D. I 
septentrionalis) . 

Materials and methods 

Taxon sampling and laboratory procedures 

A total of 26 populations from Piedmont and Liguria (Maritime Alps, and western 
Liguria Apennines in northwest Italy), were analyzed in this study (Fig. 1, Table 1). 

DNA was isolated from leg muscle of each individual, using a GenElute Mam- 
malian Genomic DNA minipreparation Kit (Sigma-Aldrich, St. Louis, MO, USA) 
resuspended in 200 |il of sterile water and stored at -40 °C. 

The entire Cytochrome Oxidase I gene (COI, total of 1500 bp), a 550-bp fragment 
of the 16S rRNA gene, a 450-bp fragment of the 12S rRNA gene and were amplified 
through the polymerase chain reaction (PGR) and sequenced from each individual. 
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Table I . North Western Italian Dolichopoda sampled in this study. All populations were analyzed at molecu- 
lar level, while only the populations in which adults were present, were analyzed also at morphological level. 



Taxon 


Localities and sample size 


Code 


0 


GenBank Accession No. 


D. azami 


Mine close to Chauves-souris Cave, 
Chateaudouble, Var, France, type locality, 
5i5 m asi; Sampled mdividuais: 3 


CHA 


X 


12S: KM086460; 16S: KM086485; 
COI: KM086535; 28S: KM086510 


D.l. ligustica 


laragnina Cave, Balestrino, oV, Liguria, close 
to Santa Lucia Inferiore Cave, 310 m asl; 
Sampled individuals: 5 


ART 




12S: KM086464; 16S: KM086489; 
COI: KM086539; 28S: KM086514 


Santa Lucia Inferiore Cave, Toirano, SV, 
Liguria, type locality, 194 m asl 


TOR 


X 




D. I. septentrionalis 


Borna Maggiore del Pugnetto Cave, 
Mezzenile, Lanzo Valley, TO, Piedmont, type 
locality, 820 m asl; Sampled individuals : 8 




V 
A 


12S: KM086474; 16S: KM086499; 
COI: KM086549; 28S: KM086524 


Dolichopoda sp. 


Tana del Peccetto Cave, Magliolo, SV 


PEC 


X 




Besta Cave, Vievola, lenda, r ranee, 910 m 
asl; Sampled individuals: 2 


BES 


X 


12S: KM086459; 16S: KM086484; 
COI: KM086534; 28S: KM086509 


Falconiere Cave, Garessio,CN, Piedmont, 
1000 m asl 


GAR 


X 




Orso Cave, Ponte di Nava, Ormea, CN, 
Piedmont, 808 m asl; Sampled individuals: 3 


GOR 




12S: KM086472; 16S: KM086497; 
COI: KM086547; 28S: KM086522 


Viozene, Ormea, upper Tanaro Valley, CN, 
Piedmont (epigean), 1250 m asl; Sampled 
individuals: 2 


LIG 




12S: KM086470; 16S: KM086495; 
COI: KM086545; 28S: KM086520 


Vernante hypogeum, Vernante, Vermenagna 
Valley, CN, Piedmont, 800 m asl; Sampled 
individuals: 5 


SVE 


X 


12S: KM086458; 16S: KM086483; 
COI: KM086533; 28S: KM086508 


Bossea Cave, Frabosa Soprana, Corsaglia 

■\ r 11 T T^' 1 O O 1 O 11 

Valley, CN, Piedmont, 836 m asl; sampled 
individuals: 2 


BOS 


X 


12S: KM086465; 16S: KM086490; 
COI: KM086540; 28S: KM086515 


Arma de Faie Cave, Albissola, Sansobbia 

A r 11 OA T T ' ' / 1 

Valley, 5V, Liguria, 624 m asl 


ARE 


X 




Bandito Cave, Valdieri, Gesso Valley, CN, 
Piedmont, 726 m asl; Sampled individuals: 2 


BAN 


X 


12S: KM086471; 16S: KM086496; 
COI: KM086546; 28S: KM086521 


Bandito hypogeum, Roaschia, Gesso 
Valley, CN, Piedmont, 730 m asl; Sampled 
individuals: 7 


SBA 


X 


12S: KM086463; 16S: KM086488; 
COI: KM086538; 28S: KM086513 


Dronera Cave, Vicoforte Mondovi, Ermetta 
Valley, CN, Piedmont, 525 m asl; Sampled 
individuals: 6 


TDR 


X 


12S: KM086468; 16S: KM086493; 
COI: KM086543; 28S: KM086518 


Rittana Cave, Rittana, Stura di Demonte 
Valley, CN, Piedmont, 1000 m asl; Sampled 
individuals: 3 


RIT 


X 


12S: KM086467; 16S: KM086492; 
COI: KM086542; 28S: KM086517 


Vallone Saben eastern hypogeum, Valdieri, 
Gesso Valley, CN, Piedmont, 750 m asl; 
Sampled individuals: 6 


svs 


X 


12S: KM086462; 16S: KM086487; 
COI: KM086537; 28S: KM086512 


Gaiola Cave, Gaiola, Stura di Demonte 
Valley, CN, Piedmont, 1020 m asl; Sampled 
individuals: 3 


PGA 


X 


12S: KM086457; 16S: KM086482; 
COI: KM086532; 28S: KM086507 
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Farout Cave, Pradleves, Grana Valley, CN, 
Piedmont, 1050 m asl; Sampled individuals: 3 


FAR 


X 


12S: KM086461; 16S: KM086486; 
COI: KM086536; 28S: KM086511 


Balmarossa Cave, Pradleves, Grana Valley, 
CN, Piedmont, 1180 m asl; Sampled 
individuals: 2 


BAL 


X 


12S: KM086452; 16S: KM086477; 
COI: KM086527; 28S: KM086502 


Buco del Partigiano Cave, Roccabruna, Maira 
Valley, CN, Piedmont, 1 100 m asl; Sampled 
individuals: 7 


BPR 




12S: KM086455; 16S: KM086480; 
COI: KM086530; 28S: 
KM086505; 


Partigiani Cave, Rossana, Varaita Valley, CN, 
Piedmont, 615 m asl; Sampled individuals: 8 


GPR 


X 


12S: KM086469; 16S: KM086494; 
COI: KM086544; 28S: KM086519 


Buco delle Locuste Cave, Pagliano, Varaita 
Valley, CN, Piedmont, 595 m asl; Sampled 
individuals: 2 


BLO 


X 


12S: KM086466; 16S: KM086491; 
COI: KM086541; 28S: KM086516 


Buco del Maestro Cave, Paesana, Po Valley, 
CN, Piedmont, 750 m asl; Sampled 
individuals: 6 


BMA 




12S: KM086453; 16S: KM086478; 
COI: KM086528; 28S: 
KM086503; 


Chiabrano Cave, Perrano, Germanasca 
Valley, TO, Piedmont, 1080 m asl; Sampled 
individuals: 6 


CHI 


X 


12S: KM086456; 16S: KM086481; 
COI: KM086531; 28S: KM086506 


Villar Focchiardo Cave, Villar Focchiardo, 
Susa Valley, TO, Piedmont, 460 m asl; 
Sampled individuals: 6 


CVF 




12S: KM086473; 16S: KM086498; 
COI: KM086548; 28S: KM086523 


Val della Torre Mine, Val della Torre, 
Casternone Valley, TO, Piedmont, 840 m asl; 
Sampled individuals: 5 


VDT2 




12S: KM086475; 16S: KM086500; 
COI: KM086550; 28S: KM086525 


Wolf Cave, Mezzenile, Lanzo Valley TO, 
Piedmont, 813 m asl; Sampled individuals: 4 


PUGl 


X 


Same as PUG 


Borna del Servais Cave, Ala di Stura, Ala 
Valley, TO, Piedmont, 1440 m asl; Sampled 
individuals: 2 


BCS 




12S: KM086454; 16S: KM086479; 
COI: KM086529; 28S: 
KM086504; 


Lacca Selva Cave, Zandobbio, BG, 
Lombardy, 340 m asl; Sampled individuals: 2 


ZAN 


X 


12S: KM086476; 16S: KM086501; 
COI: KM086551; 28S: KM086526 



The large subunit of the nuclear ribosomal DNA (28S rRNA) was also included. The 
primers used were: LCO1490, HCO 2198 (Folmer et al. 1994), UEAl, UEA5 and 
UEAIO (Lunt et al. 1996) for the COI gene, 12Sai, 12Sbi (Kocher et al. 1989; Simon 
et al. 1994) for the 12S rRNA gene and 16Sar, 16Sbr (Simon et al. 1994) for the 
16S gene. 28S rRNA was partially amplified and sequenced for a fragment of 580 
base pairs, belonging to domains 3-5, using primers from Friedrich and Tautz (1997). 
Double-stranded amplifications were performed with a Perkin-Elmer-Cetus thermal 
cycler in 25 |^1 reaction volume containing genomic DNA (10—100 ng), 1.5 mM 
MgCl2, 0.2 mM of each dNTP, 0.1 |iM primer, 1.5 units EuroTaq (Euroclone, UK) 
and the buffer supplied by the manufacturer. Optimal cycling parameters varied for 
each primer pair used. PGR products were purified using the ExoSAP digestion (Am- 
ersham Pharmacia Biotech), directly sequenced in both directions using the BigDye 
terminator ready- reaction kit, and resolved on ABI 3100 Genetic Analyzer (PE Ap- 
plied Biosystems), following the manufacturer's protocols. Sequence data were edited 
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and compiled using CodonCode Aligner (version 3.7.1). All sequences were submitted 
to GenBank (Accession Numbers are reported in Table 1). 

Each gene fragment (12S, 16S, COI and 28S) was considered separately for the align- 
ment. Sequences of 16S, 12S and 28S were aligned using CLUSTAL_X 1.81 (Thompson 
et al. 1997) with opening gap =10 and extending gap = 0.10. Cytochrome oxidase I nucle- 
otide sequences were assembled, aligned, and translated with Codon-Code Aligner 3.7.1. 



Data analysis 

Genetic differentiation and structure in the westernmost Dolichopoda species 

Haplotype and nucleotide diversity were calculated using DNASP 5.10 (Rozas et al. 
2003). Genetic distances between haplo types and between groups of haplotypes were 
calculated using p-distance as implemented in MEGA 5.2 (Tamura et al. 2007). NET- 
WORK 4.502 (Bandelt et al. 1999) was employed to calculate a median joining net- 
work representing the genealogical relationships among mtDNA haplotypes. 

A Mantel test (Mantel 1967) with 5,000 simulations was used to test for an iso- 
lation-by-distance (IBD) signature (a positive correlation between geographic and ge- 
netic distances; Wright 1943; Slatkin 1993). 

The hierarchical distribution of genetic variation was characterized using analysis 
of molecular variance (AM OVA). This method apportions genetic variation within 
and among groups, estimating F-statistics (Weir and Cockerham 1984; ExcofRer et al. 
1992; Weir 1996) that are analogous to Wright's hierarchical fixation indices under 
the island model of gene flow (Wright 1951). Three-level AMOVA was conducted 
in ARLEQUIN 3.5.1.2 (Excoffier et al. 1992; Excoffier et al. 2005) using an T^^^-like 
estimator. Samples were partitioned by geographic regions, populations within geo- 
graphic regions, and populations. These tests included permutation of inferred haplo- 
types among groups {F^^; individual haplotypes among populations but within group 
{F^^; inferred haplotypes among populations {F^^. 



Species delimitation and phylogenetic analysis 

Species delimitation analysis was carried out with two methods, a) the "classical" DNA 
barcoding gap analysis (Hebert et al. 2003) and b) the Automatic Barcode Gap Discov- 
ery (ABGD, Puillandre et al. 2012). Both analyses were carried on the present samples 
and all Dolichopoda species previously analyzed (Table 2 and Allegrucci et al. 2005, 
2011; Martinsen et al. 2009), considering only COI data set and only the common 
base pairs consisting of 964 bp.The "classical" DNA barcoding gap analysis was based 
on the pairwise genetic distances (p-distance). Distribution of genetic distance values 
was analyzed at different taxonomic levels. 
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Table 2. Species, populations, localities and accession numbers to GenBank for the Dolichopoda sequenc- 
es included in this study and retrieved from GenBank for phylogenetic reconstruction of the western Med- 
iterranean Dolichopoda species (Data from Allegrucci et al. 2005, 2009, 201 1 and Martinsen et al. 2009). 



Taxon 


Localities 


GenBank Accession No. 


Outgroup 


Euhadenoecus insolitus 


Indian Grave Point Cave (IND), De 


12S: EF216947; 16S: AY793563; 


Kalb Co., TN, USA 


COI: AY793591; 28S: EF217005 


Hadenoecus cumber- 


Bat Cave (BAT), Carter Cave State 


12S: EF216948; 16S: AY793562; 


landicus 


Park, Carter Co., KY, USA 


COI: AY793592; 28S: EF217004 


Troglophilus cavicola 


Covoli di Velo Cave (TRO), Veneto, 
Northern Italy 


12S: EF216946; 16S: AY793564; 
COI: AY793624; 28S: EF217003 


Ingroup 


D. holivari 


Forat negre Cave, Serra del Llerida, 


16S: AY507579-80, COI: AY507648-49 


Pyrenees, Spain 




D. linden 


Sirach Cave (SIR) Eastern Pyrenees, 


12S: JF826039; 16S: AY793567; 


Western-South France 


COI: AY793598-99; 28S: JF826069 




Poscola Cave (PSC), Veneto, North- 


12S: JF826054; 16S: AY793581; 




ern-East Italy 


COI: AY793611 and 13; 28S: IF826076 


D. laetitiae 


Plane Cave (GDP), Umbria, Central 
Italy 


12S: JF826053; 16S: AY793582; 

COI: AY793610 and 12; 28S: ]F826075 




Diavolo Cave (DIA), Tuscany, Cen- 


12S: JF826052; 16S: AY793580; 




tral Italy 


COI: AY793614-15; 28S: IF826074 




Valmarino Cave (VAL), Latium, 


12S: JF826055; 16S: AY793583; 




Central-South Italy 


COI: AY793616-17; 28S: IF826077 




Fontanelle Cave (FON), Campania, 


12S: JF826057; 1 6S: AY7935 84; 


D. geniculata 


Southern-West Italy 


COI: AY793594; 28S: IF826079 


Ischia cellars (ISC), Ischia Island, 


12S: JF826058; 16S: AY793585; 




Campania, Southern-West Italy 


COI: AY793595; 28S: IF826080 




Roman Aqueduct (PNZ), Ponza 


16S: AY793586; COI: AY793596-97 




Island, Latium, Southern-West Italy 




D. capreensis 


San Michele Cave (CPR), Capri Is- 


12S: JF826059; 16S: AY793587; 


land, Campania, Southern Italy 


COI: AY793606-7; 28S: IF826081 


D. palpata 


Tremusa cave (TRE), Calabria, 


12S: JF826060; 16S: AY793588; 


Southern Italy 


COI: AY793608-9; 28S: IF826082 


D. baccettii 


Punta degli Stretti Cave (PST), Tus- 


12S: JF826046; 16S: AY793571; 


cany, Central-West Italy, 


COI: AY793639-40; 28S: IF826068 


D. aegilion 


Campese Mine (CAM), Giglio Is- 


12S: JF826045; 16S: AY793570; 


land, Tuscany, Central-West Italy 


COI: AY793600; 28S: IF826067 




Pipistrelli Cave (ORS), Tuscany, 


12S: JF826044; 16S: AY793573; 




Central-West Italy 


COI: AY793633; 28S: IF826066 


D. schiavazzii 


Marciana Cave (MRC), Elba Island, 


12S: JF826043; 16S: AY793572; 


Tuscany, Central-West Italy 


COI: AY793635; 28S: JF826065 




Fichino Cave (FIG), Tuscany, Cen- 


12S: JF826042; 16S: AY793574; 




tral-West Italy 


COI: AY793634-36; 28S: IF826064 




Brando Cave (BRA), Corsica Island, 


12S: JF826047; 1 6S: AY793578; 




France 


COI: AY793627 and 31-32; 28S: 


D. bormansi 




IF826069 


Sisco Cave (SIS), Corsica Island, 


12S: JF826048; 1 6S: AY793579; 




France 


COI: AY793625-26 and 28; 28S: 
IF826070 
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D. cyrnensis 


Valletto Cave (VLl ), Corsica Island, 
France 


12S: JF826050; 16S: AY793577; 
COI: AY793620-21; 28S: IF826072 


Sahara Cave (SAB), Corsica Island, 
France 


12S: JF826049; 16S: AY793576; 
COI: AY793618-19; 28S: IF826071 


D. muceddai 


Limbara Mount, Sardinia Island, Italy 


12S: JF826051; 16S: AY793575; 
COI: AY793629-30; 28S: IF826073 



The ABGD method (Puillandre et al. 2012) automatically finds the distance at 
which a barcode gap occurs and sorts the sequences into putative species based on this 
distance. The method statistically infers the barcode gap from the data and partitions 
the data accordingly. Populations belonging to the same species therefore should be 
grouped in the same partition. This procedure is then recursively applied to the previ- 
ously obtained groups of sequences. COI alignments were uploaded at http://wwwabi. 
snv.jussieu.fr/public/abgd/abgdweb.html and ABGD was run with the default settings 
(Pmin = 0.001, Pmax = 0.1, Steps = 10, X (relative gap width) = 1.5, Nb bins = 20) 
and using p-distance. 

A phylogenetic analysis was carried out considering the four sequenced genes, by 
constructing a concatenated matrix, partitioned by genes. Phylogenetic analysis was 
performed only on the western Mediterranean Dolichopoda species, using Bayesian in- 
ferences as implemented by the software MrBayes 3.1b4 (Huelsenbeck and Ronquist 
2001). Mrmodel test (Nylander 2004) was used to perform hierarchical likelihood 
ratio test and calculate approximate Akaike Information Criterion (AIC) values of the 
nucleotide substitution models for each gene fragment. 

At least two simultaneous searches were conducted comprising four Markov 
chains started from a randomly chosen tree and run for 1,000,000 generations, with 
sampling every 100 generations. The following descriptors were assumed to indicate 
convergence on a common phylogenetic topology by separate Bayesian searches: 
similarity in log likelihood scores at stationarity, similarity in consensus tree to- 
pologies and PP values for supported nodes, and a final average standard deviation 
of split frequencies (ASDSF) for simultaneous searches approaching zero. The first 
1,000 trees were discarded as burn-in and posterior probabilities (PP) were calcu- 
lated from post burn-in trees. 



Morphological analysis 

The most important morphological characters used in the systematics of the ge- 
nus Dolichopoda are the shape of the epiphallus and of the male tergum 10. In 
this study, the investigation of these characters was integrated by the analysis of 
other morphological features as male sub-genital plate, female sub-genital plate and 
shape of ovipositor. 

Overall, we have checked 22 populations from the northwestern Italy for a total of 
64 adult individuals. All these populations are the same analyzed for the genetic varia- 
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tion (Table 1). Further comparisons were made using three more samples (PEC, ARF, 
GAR) coming from the collection of Museum of Natural History of Genoa G. Doria. 
The specimens were studied using a stereomicroscope Leica MZ 12.5. All the measures 
of the morphological parameters are in mm. Photos were taken with a digital camera, 
Nikon Coolpix P300. Pictures were processed using a digitiser board WACOM CTH 
461 and Adobe Photoshop CS3 Extended Version 10.0. 

Table 3 shows the measures of twelve morphological characters in 6 populations 
where a sufficient number of adults were available. In particular, we considered the to- 
tal body length, the pronotum length, and the length of tarsum, femur and tibia of the 
fore, middle and hind legs in both males and females and the ovipositor length in fe- 
males. Measurements were taken using a digital caliper (accuracy 0.01 mm). Morpho- 
metric variables were analyzed by using the non-parametric ANOVA Kruskal-Wallis 
test. A multiple regression analysis was carried out comparing morphometric variables 
with the altitude of the caves. 



Results 

Genetic difFerentiation and structure in the westernmost Dolichopoda species 

A total of 2450 mitochondrial base pairs corresponding to the entire CO/ gene, to 450 
bp of 12S and 550 bp of 16S were successfully sequenced and aligned in 104 specimens 
of Dolichopoda populations from northwestern Italy. 28S rRNA, consisting of 500 bp 
was also sequenced in all considered samples. Ribosomal genes were identical in all as- 
sayed samples. Therefore, the analysis was limited to the COI mitochondrial gene. 

The 104 individual COI sequences comprised a total of 15 haplotypes determined 
by 13 segregating sites. Global haplotype diversity, h, was 0.875 and nucleotide diversity, 
/>, was 0.00231. 

Dolichopoda haplotypes were organized in three main haplogroups in which haplo- 
types are ordered geographically (Fig. 2) . In particular, haplogroup 1 included samples 
from the northernmost area, haplogroup 2 samples from the westernmost area and 
haplogroup 3 consisted of samples from the south-eastern area (Fig.2). 

To investigate this further we carried out a genetic distance analysis using COI 
gene as a barcode and p-distance between all studied Dolichopoda species (Allegrucci 
et al. 2005, 2009, 201 1; Martinsen et al. 2009). Genetic distance values at intra- and 
interspecific levels are compared in Figure 3. In particular. Figure 3 shows the genetic 
distance values found in intra- and interspecific comparisons of all studied Dolichop- 
oda species compared to the interspecific values found in comparisons between the 
northwestern Italian populations. Intraspecific values between all Dolichopoda species 
ranged from 0 to 0.023, with a mean of 0.001 (± 0.0019, standard deviation), while in- 
terspecific values ranged from 0.016 to 0.134, with a mean of 0.134 (± 0.020, standard 
deviation). Pairwise comparisons between the northwestern Italian populations ranged 
from 0 to 0.007 with a mean of 0.003 (± 0.0015, standard deviation). 
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Figure 2. Median-joining network analysis in the populations of Dolichopoda considered in this study: 
circled areas are proportional to the number of individuals sharing the same haplotype; bars, along con- 
nections, indicate the number of nucleotide substitutions. Individuals from the same population can 
show different haplotypes. Based on number of nucleotide substitutions, three haplogroups can be distin- 
guished, each indicated with different color. 



ABGD analysis proposed several partitions that varied according to the different a 
priori thresholds. Apart from the two extreme a priori threshold values {P = 0.001 and P 
= 0.035), for which aberrant number of species hypotheses were obtained (almost every 
haplotype was considered as a different species hypothesis for P = 0.001 and, conversely, 
all the haplotypes were combined in a single species hypothesis (or P = 0.035), all the 
tested a priori thresholds lead to the same splitting. The westernmost Dolichopoda taxa 
were grouped together, while the other groups corresponded to the nominal species. 

Analysis of molecular variance (AMOVA; Excoffier et al. 1992) used to measure 
the proportion of genetic variation among subdivided populations, suggested a strong 
structure within the populations {F^^ = 0.979, P=0). Genetic variation among different 
geographic groups and among populations within each group was 52.16 % and 45.8 
%, respectively, with F^^ = 0.522 {P = 0) and F^^ = 0.957 {P = 0). Mantel test sug- 
gested a clear isolation by distance across the sampled region (R^ = 0.623, P < 0.002). 
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Figure 3. Distribution of genetic distance values (p-distances) at different taxonomic levels. Pairwise 
comparisons at the intra- and inter-specific level in all Dolichopoda species are reported. 



Phylogenetic analysis 



MRMODELTEST (Nylander 2004) indicated GTR + I + G (Lanave et al. 1984; Gu 
et al. 1995) and K80 + I (Kimura 1980) as the best models of DNA substitution for 
the mitochondrial genes and the nuclear one, respectively. The phylogeny in Figure 4 
supports the major phylogenetic relationships previously demonstrated (Allegrucci et 
al. 2005, 201 1; Martinsen et al. 2009). The Italian Dolichopoda species are separated in 
two main clusters, as expected. The first group included the continental species from 
Liguria to southern Italy {D. ligustica, D. I. septentrionalis, D. azami, D. laetitiae and 
D. geniculata) . Populations from D. ligustica, D. I. septentrionalis, and D. azami consti- 
tuted a homogeneous clade well differentiated from D. laetitiae and D. geniculata. The 
second group was represented by the southernmost species from Tuscany, to Corsica, 
Campania and Calabria (i.e., D. aegilion, D. baccettii, D. schiavazzii, D. cyrnensis, D. 
bormansi and D. muceddai and D. palpata with its sister taxon D. capraensis) . Linked 
to these clusters are D. linderi and D. bolivari. 
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Figure 4. Relationships among Italian species of Dolichopoda inferred from Bayesian analysis based on 
three mitochondrial and one nuclear genes. Values above branches indicate posterior probabilities derived 
from Bayesian analysis. Scale bars: 0.03 substitutions per site. Only posterior probability (PP) values > 
0.80 are shown. 
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Figures 5-13. Morphological characters of Dolichopoda azami (CHA type locality), Dolichopoda I li- 
gustica (TOR type locality) and Dolichopoda I. septentrionalis (PUG type locality). Left: dorsal view of 
tergum 10; centre: epiphallus of the male; right: particular of the female sub-genital plate. Arrows show 
diagnostic characters. 

Morphological analysis 

In Figures 5-13 morphological comparisons among the three samples from the typical 
locality of D. azami (Chateaudouble, CHA), D. ligustica (Toirano, TOR) and D. I. sep- 
tentrionalis (Pugnetto, PUG) are reported. In these figures, the three main characters are 
compared: male tergum 10, epiphallus and female sub-genital plate. The other two con- 
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Figures 1 4-22. Morphological characters of the Dolichopoda from Liguria: ARF, GAR and PEC. Left: 
dorsal view of tergum 10; centre: epiphallus of the male; right: particular of the female sub-genital plate. 
Arrows show diagnostic characters. 

sidered morphological features, male sub-genital plate and ovipositor, are not reported 
because showed a very uniform pattern. As already evidenced in the taxonomic revision 
of the genus Dolichopoda (Baccetti and Capra 1959), it is possible to distinguish the 
species D. azami through the male tergum 1 0 that shows two rounded lobes separated 
by a large concavity with two little acuminate tips (Fig. 5) and median process of the 
epiphallus enlarged at the basis (Fig. 6). In the female, the sub-genital plate is completely 
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rounded with a narrow apical incision (Fig. 7). On the contrary, D. ligustica shows a 
male tergum 10 with little oval lobes separated by almost straight margin (Fig. 8), and 
median process of epiphallus pyramidal with lateral margin almost straight (Fig. 9). The 
female sub-genital plate is very similar to that of D. azami (Fig. 10) but the median 
incision appears less pronounced. Finally the subspecies D. I. septentrionalis is distin- 
guishable only for the more slender median process of the epiphallus (Fig. 12). Based on 
these observations, it is evident the great similarity of the three taxa here considered and 
consequently the strong difficulty in their morphological distinction. 

In Figures 14—22, we compared three samples from western Liguria (Table 1). The 
three considered characters appear uniform resembling all to the morphology of TOR 
(Figs 8-10) the typical D. ligustica. On the contrary, the analysis of a series of samples from 
northern and southern Piedmont (Figs 23-37) highlights a more complex situation with 
the contemporary presence of characters typical of one or the other species in the same indi- 
vidual. In particular, even if all the specimens analyzed show the tergum 10 very similar to 
that of D. azami (cfr Fig. 5 with Figs 23, 26, 29, 32, 35) the same specimens often show the 
shape of epiphallus and of the female sub-genital plate more similar to that of D. ligustica. 

In Table 3, the measures of 12 morphological characters, recorded on 6 popula- 
tions, where a sufficient number of adults were available, are reported. Mann-Whitney 
test carried out between males and females morphometric traits revealed significant 
differences {P < 0.05), being the females smaller than males. ANOVA analysis high- 
lighted significant differences between some morphometric measures in specimens 
coming from caves situated at different altitude both in males and in females. Mul- 
tiple regression analysis revealed a significant correlation between the morphometric 
measures and the altitude of the cave (R^ = 0.661, Wilk's Ibd = 0.0045; P « 0.01). In 
particular, samples coming from caves at low altitudes resulted bigger than those from 
caves at higher altitudes. 



Discussion 

The Northwestern Italian Dolichopoda species group comprises a uniform assemblage 
of populations, both genetically and morphologically. From genetically point of view, 
the median joining network analysis (Fig. 2) evidenced three main haplogroups geo- 
graphically distributed, but the resulting scenario is rather homogeneous with a maxi- 
mum number of mutation events equal to 8. Genetic distance values based on COI 
gene revealed a very little differentiation between all studied samples. All pairwise 
comparisons fall in the range of intraspecific comparisons observed in all the other 
Dolichopoda species, with a mean of 0.3% (Fig. 3). These results suggest a level of dif- 
ferentiation usually found between populations belonging to the same species rather 
than to different species. ABGD analysis confirms these results, grouping together all 
the northwestern Dolichopoda taxa. 

However, despite the low variability, these populations show a significant level 
of genetic structure. Indeed, AMOVA analysis evidenced significant partitioning of 
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variation within and among populations and among regions and a low level of gene 
flow between the different regions. Isolation by distance was also detected by Mantel 
test, suggesting that dispersal rates in these Dolichopoda populations are very low 
and limited to the geographically close ones. This result is not surprising for cave 
dwelling organisms because it reflects the geographic isolation of the different caves 
and / or groups of caves to which Dolichopoda populations are confined (Sbordoni 
et al. 2012). The Maritime Alps region is characterized by a combination of Medi- 
terranean, Apennine and Alpine flora and this diversity in habitat geography could 
strongly influence the different degree of isolation of the Dolichopoda populations. 
While the presence of beech and oak forests would account for continuing gene flow 
among Dolichopoda species, the occurrence of Alpine forests and to some extent 
of the Mediterranean "macchia" could prevent the gene exchange (Allegrucci et al. 
1997; Sbordoni et al. 1985, 2012). 

The phylogeny illustrated in Figure 4, based on three mitochondrial (COI, 12S 
and 16S) and one nuclear (28S) genes, supports the major phylogenetic relationships 
previously demonstrated (Allegrucci et al. 2005, 2011) highlighting that the north- 
western Italian populations represent a single, uniform, well supported cluster. With 
respect to the other Italian species, they are genetically very similar, showing very short 
branch lengths and forming a single cluster in which the relationships between the 
different populations are not clearly resolved. Based on previous studies carried out 
on the 90% of the species belonging to genus Dolichopoda (Allegrucci et al. 2005, 
201 1), we know that divergence time between these populations and their sister taxon, 
represented by the cluster D. laetitiae-D. geniculata, is dated back to about 2 Million 
years ago, during the late Pliocene. In that time, the western Mediterranean lineages 
probably arose from a single colonization event from the North of Italy. An ancestral 
group spread in two directions, colonizing on one hand Piedmont, Liguria and France 
and, on the other side, the entire inland of Italian Peninsula. 

The morphological analysis shows high degree of similarity between all analyzed 
samples, suggesting that the morphological traits commonly used in the taxonomy of 
the genus Dolichopoda (tergum 10, epiphallus shape in the male and sub-genital plate 
in the female), are not able to discriminate clearly specimens belonging to D. azami, 
D. ligustica and D. I septentrionalis. In particular, although very few morphological 
differences can be observed between specimens from the typical localities of the three 
taxa (CHA, typ. loc. of D. azami\ TOR, typ. loc. of D. ligustica and PUG, typ. loc. 
of D. I. septentrionalis), the most of the studied samples show some variability in the 
diagnostic characters occurring more or less randomly. Often we observed in the same 
specimen, the contemporary presence of characters typical of one or the other taxon. 

Morphometric analysis suggested significant differences both between males and 
females and between samples from caves situated at different altitudes. Caves situ- 
ated at high altitude are presumably subjected at lower temperatures and therefore 
the bioclimatic factors may be the major determinants of the morphometric patterns 
observed. Morphological traits might be subjected to substantial non-genetic variation 
and/or might show local short-term adaptations to peculiar environmental conditions. 



20 



Giuliana Allegrucci et al I ZooKeys 437: 1—23 (2014) 



These differences were already been observed in populations of D. laetitiae and D. 
geniculata in the central southern Italy (Allegrucci et al. 1987; Cesaroni et al. 1997). 
In those studies, we demonstrated that variations in morphometric measures (and in 
particular in leg elongation) is almost entirely under the control of an environmental 
gradient, synthetically described by the cave temperature. 

In conclusion, the scenario, depicted by both molecular markers and morphologi- 
cal traits, suggests that the northwestern Italian Dolichopoda constitute a single species 
with populations characterized from high similarity between each other and, due to 
the priority of description, they should be assigned to D. azami Saulcy, 1893 which 
is the older available name at the species rank. However, in order to portrait the geo- 
graphical partitioning and the related genetic differences we suggest to preserve the 
names ligustica Baccetti & Capra, 1959 and septentrionalis Baccetti & Capra, 1959 
to indicate groups of populations corresponding to the different COI haplogroups as 
outlined in the network analysis (Fig. 2). 

Therefore, the new taxonomic arrangement here proposed is: 

D. azami azami SmxAcyy 1893 

Type locality. Mine close to Chauves-souris Cave, Chateaudouble, Var, France, 525 m asl; 
D. azami ligustica Baccetti & Capra, 1959, stat. n. 

Type locality. Santa Lucia Inferiore Cave, Toirano, SV, Liguria, 194 m asl; 
D. azami septentrionalis Baccetti & Capra, 1959, stat. n. 

Type locality. Borna Maggiore del Pugnetto Cave, Mezzenile, Lanzo Valley, TO, 
Piedmont, 820 m asl. 
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